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Magnetars^ comprise two classes of rotating neutron stars (Soft Gamma Re¬ 
peaters (SGRs) and Anomalous X-ray Pulsars), whose X-ray emission is powered 
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by an ultrastrong magnetic field, B ~ 10^® G. Occasionally SGRs enter into active 
episodes producing many short X-ray bursts; extremely rarely (about once per 
50 years per source), SGRs emit a giant fiare, an event with total energy at least 
10^ higher than their typical bursts.^ ^ Here we report that, on 2004 December 
27, SGR 1806—20 emitted the brightest extra-solar transient event ever recorded, 
even surpassing the full moon brightness for 0.2 seconds. The total (isotropic) 
fiare energy is 2 x 10^® erg, ~100 times higher than the only two previous events, 
making this fiare a once in a century event. This colossal energy release likely 
occurred during a catastrophic reconfiguration of the magnetar’s magnetic field. 
Such an event would have resembled a short, hard Gamma Ray Burst (GRB) 
if it had occurred within 40 Mpc, suggesting that extragalactic SGR fiares may 
indeed form a subclass of GRBs. 

Only two other giant flares have previously been recorded, one each from SGR 0526-66 
on 1979 March 5^’^ and SGR 1900-1-14 on 1998 August 27.^ Intense X-ray burst activity 
from SGR 1900-1-14 preceded the 1998 August 27 flare^; there are very few observations of 
SGR 0526-66 preceding the 1979 March 5 event. In the year leading up to the SGR 1806-20 
flare, well-sampled X-ray monitoring observations of the source with the Rossi X-ray Timing 
Explorer (RXTE) indicated that it was also entering a very active phase,® emitting more 
frequent and intense bursts and showing enhanced persistent X-ray emission that was, indeed, 
a prelude to the unprecedented giant flare. 

On 2004 December 27 the Swift satellite^ was among a large number of spacecraft inun¬ 
dated by radiation from SGR 1806-20.®’®’®® The Burst Alert Telescope (BAT)®"^ is a 7-ray 
(15-350 keV) coded aperture imager on Swift. Although Swift was turned away from the SGR 
location, and so the event illuminated the detector from behind, the flux that passed through 
the spacecraft and shielding of the BAT provided excellent measurements of the event. The 
BAT light curve (Figure 1, and with more detail in the supplementary figures) demonstrates 
that magnetar giant flares are remarkably similar: all three start with an initial very short 
and spectrally hard main spike, followed by an extended softer tail highly modulated at the 
neutron star’s spin period. The bright, main spike lasts ~ 0.5 s and is followed by a tail 
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with ~ 50 cycles of high-amplitude pulsations at the known rotation period of SGR 1806-20 
(7.56s). In the Dec. 27 event we also notice a 1-s long, flat-topped precursor burst at 142 s 
before the main spike. 

Several astounding new properties of a magnetar flare are revealed from the superb time 
resolution of the BAT. Figure lb plots the sharp initial rise of the main spike in time bins 
of lOOps, equivalent to the light crossing time of the neutron star diameter. Prior to the 
steep rise of the initial spike, the count rate was rising for 40 ms at a slower rate (shown in 
supplementary figure 2) and had reached roughly 30,000 cps (above a ~ Okcps background) 
by t = 0. At that point it increased by a factor of more than 100 in less than 1.5 ms 
(0.3ms e-folding rise time). This is followed by at least one dip and continued brightening 
(additional dips would not be visible due to instrument saturation) on its way to the peak. 
The flare rise is thus resolved for the first time. The flux during the spike, though heavily 
attenuated, saturated the BAT modules, precluding a reliable flux measurement. We have 
used, therefore, the SOPA®^ and ESP®^ instruments located on geosynchronous satellites (see 
Supplementary Methods Section) to measure the main peak flux. The SOPA instruments are 
small silicon detectors designed with fast event processing to measure the high particle fluxes 
found in orbit. During the peak of the burst, each detector had a dead-time greater than 
50%, but this level of saturation can be accurately corrected for. We £t the SOPA data with 
an exponential-cutoff power law (finding a characteristic temperature kT = 0.48(4) MeV and 
a power-law photon index —0.2(1)) and derive a flux of 5.0(3) erg cm“^ s“^ over an 0.160 s 
integration time for 45keV to 10 MeV photons; the corresponding fluence is 0.80(5) ergcm”^. 
This duration and spectral hardness is in the range of characteristics found for the short, 
hard subclass of classic GRBs.^^ 

Since the count rate was significantly lower during the tail, we were able to model the off- 
axis illumination and calibrate the flux and spectroscopy measurements (see Supplementary 
Methods). We find the tail of the burst to have an energy fluence of 1.0(5) x 10“^ ergcm“^ 
at photon energies > 60 keV. The spectral fits are consistent with thermalized spectra with 
kT ~ 15 — 30 keV as seen in previous flares, implying a comparable energy fluence below our 
60keV threshold. Accounting for the 10 — 60 keV photons, we project the total tail fluence to 
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be ^ 2 X 10“^ ergcin”^, roughly 0.3% that of the main peak. For a distance to SGR 1806-20 
of ISdiskpc/® we then hnd an (isotropic equivalent) energy release of 2 x lO^^d^gerg in the 
spike, and 5 x erg in the tail. Thus, the isotropic equivalent energy in the initial 

spike is about two orders of magnitude larger than that in the other two giant flares, while 
the energy in the tail is comparable. Indeed, a radio afterglow^® was detected from this 
flare with a luminosity 500 times that from the August 27 flare, suggesting a truly large 
difference in the prompt burst energy. The consistency of the tail energies among the three 
flares is attributable to the fact that they are limited by the storage capacity of the magnetic 
field^’^^ and should be as constant from source to source as the held energy. Thus, the tail 
luminosities, which are expected to be at roughly the magnetically modihed Eddington limit, 
should also be similar, as observed. The extent of magnetic reconnection, on the other hand, 
governs the prompt energy release during the main spike; this can vary greatly from one 
event to the next, even within the same source. 

The pulse prohle in the tail of the hare just after the main spike features one large 
peak and two smaller adjacent local maxima separated by about 1/4 of a rotation cycle 
(Figure 2 and supplementary figure 1). The relative intensities of the peaks change during 
the tail, but their phases remain hxed, indicating that the held conhguration does not change 
substantially during the tail and that the released energy comes from the trapped hreball. 

The polar R—held of SGR 1806-20 has been calculated^® from its spin-down rate to be 
~ 1.6 X 10^^ G, corresponding to a external magnetic held energy of 2 x 10^^ erg, which 
indicates that at most lOdi^/ f such giant hares can be produced from the star in its lifetime 
(here / is the beaming factor). 

We used RXTE to measure the spin frequency and spin-down rate of the SGR 30 days 
after the hare.^® The frequency is consistent with an extrapolation of the pre-hare frequency 
with pre- and post-hare spin-down rates. Thus, the December 27 hare could not have caused a 
rapid, lasting change in the spin frequency greater than ~ 2 x 10“® Hz; this, despite the much 
larger apparent burst energy, limits the frequency change to be at most comparable to that 
seen following the August 27 hare.^'^ The post-hare spin-down rate, —3.15(9) x 
although lower than it was shortly before the hare, is still in its historical range. 
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The three time scales in the phenomenon, a) the rise time of ~ 1 ms, b) the duration of the 
hard spike, ~ 0.5 s, and c) the duration of the tail, several minutes - are similar for all three 
giant flares. These are attributed to the Alfven propagation times in a) the magnetosphere 
and b) the star, and c) the cooling time of the trapped pair hreball, respectively. 

Violent energy dissipation can occur anywhere in the magnetically dominated region, 
which includes the outer layers of the neutron star: if an energy of erg is dissipated 

roughly uniformly in the reconnection region of volume lO^^Vigcm^, then matter above the 
layer at density 1 x Itf'E^Q/Visg/crv? will have an energy density larger than its gravita¬ 
tional potential and become unbound. This is about 10^^ g, which can be ejected into the 
magnetosphere at fractions of the speed of light. Such a mass ejection (which need not be 
isotropic) is enough to power the observed radio nebula and its 0.3c expansion.^® 

The two earlier flares would have been detectable by existing instruments only within 
~ 8 Mpc, and, therefore, it was not previously thought that such flares could be the source of 
the short, hard GRBs. The main spike of the December 27 giant flare would have resembled 
a short, hard GRB had it occurred within 40 dis Mpc, encompassing even the Virgo cluster. 
Magnetar formation rate is expected to follow the star formation rate, which is (for z = 0) 
1.3±0.2 Mq yr“^ in our Galaxy^^ and 0.013j;o;oo7 Mq Mpc ^ yr ^ averaged over intergalactic 
scales.This suggests that the Burst And Transient Source Experiment (BATSE) onboard 
the Gompton Gamma Ray Observatory would have triggered on such events as short GRBs 
at a rate of Vbatse = 80(Ai"gai,GF/0.03 yr^^jd^s yr“^ to be compared to an estimate of the 
dvrsr BATSE rate of about 150 yr“^. Here, Agai,GF is the average rate of giant flares in the 
Galaxy similar to the Dec. 27 event. The observed isotropic distribution of short BATSE 
GRBs on the sky and the lack of excess events from the direction of the Virgo cluster suggests 
that only a small fraction (< 0.05) of these events can be SGR giant flares within < 40 Mpc. 
This implies either; d ^ 7kpc; Vgai.GF ~ (3) x 10“^ yr“^ on average for a Galaxy like our 
own; or the luminosity distribution includes even larger SGR flares that can be seen at a 
greater distance.One possible distinction of these from the classic GRB population may 
well come from their radio observations, as their radio afterglows should not be detectable 
beyond Mpc. The fraction of SGR events among what are now classihed as short GRBs 
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may not be predominant, but it should be detectable. This will be testable with future Swift 

observations. 
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Supplementary Methods 

SOPA and ESP peak fluence determination 

Each Synchronous Orbit Particle Analyzer (SOPA®^) consists of three identical particle 
telescopes that function as AE-E particle spectrometers. Energy depositions in a 3000 
fim thick, 0.25 cm^ silicon E-sensor that occur without a hnite AE signal are identihed as 
electrons and pulse-height analyzed into 10 deposited energy channels. To provide accurate 
measurements of large particle fluxes, each SOPA telescope has a small held of view (11°) and 
fast readout electronics (5.5 /xs deadtime). These characteristics also make the SOPA electron 
channels ideal for intense y-ray huxes that would overwhelm larger detectors designed for 
maximum sensitivity to weak signals. The 10 SOPA electron channels cover the deposited- 
energy range 45keV to >1.5 MeV. These data are continuously multiplexed over individual 
telescopes with a cadence and resolution of 0.160 s. The 15 SOPA telescopes on the hve 
satellites point in different directions, and the phase of the sampling time varies from satellite 
to satellite, providing diversity in observing conditions that allows an estimate of systematic 
errors by comparing measurements across instruments. 

The SOPA data for each well-placed detector were deadtime-corrected (with typical dead- 
times of ~ 50%) and ht to a spectral form consisting of a power law with an exponential 
cutoff (a form commonly used to ht gamma-ray burst and SGR spectra). For some of the 
detectors, the intense portion of the burst was split between two 160 ms time bins, in others 
the burst was almost entirely contained in a single time bin. This implies that the brightest 
part of the main spike was less than, but comparable to, 160 ms. 

The derived spectra from the spike have ‘temperature’ parameters (determining the ex¬ 
ponential cutoh energy) of kT ~ 500 keV with a spread among instruments of ~100keV, 
and power-law index parameters of —0.2 with a spread of 0.3. The exponential cutoh is 
well below SOPA’s upper energy limit. An example of such a ht is shown in Supplementary 
Figure 3. The 1-timebin or 2-timebin integrated huxes for each SOPA instrument gives 
5ergcm“^ s“^ x 0.160s, or a huence of 0.8ergcm“^ above 45keV. The spread on the derived 
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fluence parameter across SOPAs is 0.05 erg cm“^, which we take as a measure of the system¬ 
atic uncertainty. (The statistical uncertainty is miniscule compared to the systematics.) 

This flux is much higher than our experience with the previous two giant SGR flares. 
For perspective, the light from the full moon has an energy flux of 3ergcm“^ s“^. Therefore, 
verihcation from other instruments is desired to add to the measurement’s credibility. 

The Energy Spectrometer for Particles (ESP®^) is a higher-energy-range instrument on 
the same hve satellites as SOPA (a pair of ESPs flies on each satellite). Each is a scintillator- 
photomultiplier combination, with the scintillator being a 3.81 cm dia. x 3.81 cm bismuth 
germanate cylinder completely surrounded by 0.63-cm thick plastic scintillator as a phoswich. 
The operation of this instrument includes the measurement of the anode current, sampled 
every 128 ms (but with an effective integration time of a few ms). These anode currents 
have been calibrated on the ground by exposure to a ®°Co y-ray source at flux levels up to 
360ergcm“^ s“^. The highest anode currents currents measured during the superflare corre¬ 
spond to a peak flux of 6 ergcrn”^ s“^, and so are consistent with the SOPA measurements. 

Preliminary results reported®^ for the RHESSI telescope give lower limits on the fluence 
> 0.3ergcm“^ s“^ over a 3keV-15MeV energy range, and > 0.1 ergcrn”^ s“^above 65keV. 
These lower limits are consistent with our (somewhat higher) measurements. 

BAT tail fluence determination 

The Burst Alert Telescope (BAT)®"^ on Swift is a wide (1.5 sr) held of view (FOV) coded 
aperture telescope designed to rapidly identify and locate Gamma Ray Bursts (GRBs). The 
detector array is composed of 32768 individual 4x4x2 mm^ GdZnTe crystals, each with 
its own ampiher-discriminator chain, yielding a 5200 cm^ position-sensitive detector plane 
with 4 mm resolution and a nominal 15-350 keV energy range. 

As Swift was oriented at the time of the main spike, the SGR illuminated the detector at 
an angle of 105° from the center of the held of view. At this angle, the y-rays passed through 
or were scattered by the body of the spacecraft, the tube of either the X-Ray Telescope (XRT) 
or UV and Optical Telescope UVOT (depending on detector location) and a graded-Z shield 
behind the detector plane. After the first and second slew, giving angles of 61° and 67°, 
respectively, the y-rays no longer had to pass through the spacecraft body (although they 
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still encountered the XRT or UVOT and graded-Z shielding on the sides of the BAT before 
reaching the detector). 

Using an instrumental and spacecraft mass model and Monte Carlo methods, we gener¬ 
ated Detector Response Matrices (DRMs) for the 105° and 67° orientations, allowing us to 
use the XSPEC package to derive fluxes and spectra from count rates. The effective area 
of the instrument in the 105° orientation is a few hundred cm^ at photon energies above 
100 keV, dropping to ~10cm^ around 60 keV and insignihcant below 45 keV. In the 67° 
orientation the effective area was roughly twice that for 105° at photon energies of 60 keV 
and above, and maintained a small but appreciable effective area of 3 cm^ at our 15 keV 
threshold. (In that orientation, low energy photons could pass through the coded aperture 
mask and scatter off the innermost (copper) layer of the shield directly into the detector.) 
Strong instrumental fluorescence lines are predicted by the Monte Carlo, and are indeed the 
most obvious features of the measured spectrum throughout the tail emission. 

One of the data products that BAT generates in response to a GRB or other transient is 
a listing of every detected photon, including the measured energy deposition and the time to 
100/is resolution. For transients that are not located {e.g., as in the case of the SGR, when 
they are not in the BAT’s imaged £eld-of-view), 20 seconds of such data is returned. This 
is sufhcient to produce the high time resolution plot of Figure lb, but did not include most 
of the tail emission. However, continuous data is available as count rates in 4 energy bins at 
64 ms resolution. 

We integrated the count rates in the two highest energy bins (50-100 and 100-350 keV 
deposited energy) to get the background-subtracted two-channel ‘spectra’ in three time in¬ 
tervals: before, during, and after the hrst slew. 

These spectra were then fit with XSPEC using two spectral forms: an unbroken power- 
law (PL) and an Optically Thin Thermal Bremsstrahlung (OTTB). The first interval was fit 
using the 105° DRM, the third interval was £t using the 67° DRM, and the second interval, 
the slew, was fit using both. As each individual fit is constrained (two free parameters 
fitting two spectral bins), and because systematic uncertainties are large while uncertainties 
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in counting statistics are practically non-existent, fit quality is inapplicable and error 
estimates must come from arguments such as consistency levels among fits. 

The PL hts gave large spectral indices (—5-7), which is usually an indicator that 

an exponential cut-off is required by the £t. The OTTB (a spectral form that includes an 
exponential cut-off) gave fits with temperature parameters in the 15-35 keV range, consistent 
with previous measurements of giant flare tails. The fit temperature is unreliable, based solely 
on the ‘color’ (two-bin count ratio) at energies well above the peak, but the consistency gives 
qualitative support to the use of this spectral form. 

The calculated energy fluxes based on these two spectral forms were consistent within 20% 
over 60-350 keV. Fits made over the slew interval with the two DRMs varied in calculated 
fluence by a factor of 2, consistent with the change in effective area over this time interval. 
Integrating the fluences over all three intervals (using the average fluxes of the two DRMs 
during the slew interval) gave 0.93 x 10“^ and 1.1 x 10“^ ergcm“^ for the OTTB and PL, 
respectively. 

Based on the obvious systematic uncertainties in this procedure, we adopted a systematic 
error estimate of 50%. This is much larger than the discrepancy between the two models’ 
hts, but many of the systematics of fitting two-parameter models with two measurements 
are correlated. 

In addition, the width of the pulses in each cycle are narrow, indicating that the tail 
emission is substantially beamed. Measurements at Earth show 3 or 4 peaks, but observers 
at other latitudes relative to the neutron star would see a different pulse pattern, with a 
different overall intensity. {E.g., the August 27 1998 tail from SGR 1900-1-14 showed 4 main 
peaks with 8 distinct maxima per cycle.) Thus, even perfect measurements would provide 
limited advantage in determining the total energy radiated during the tail. 

Averaging the fluences derived from the two spectral forms yielded our stated fluence of 
1.0(5) X 10“^ergcm“^. 


BAT deadtime correction 
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bat’s 32,768 individual CdZnTe detectors are grouped into 256 units, called ’sand¬ 
wiches’. Each sandwich contains 128 detectors feeding an integrated circuit (the XA-1) that 
contains 128 ampliher-discriminator chains and a demultiplexer. When a detector absorbs 
a 7 -ray, the signal is amplihed, then detected by the discriminator, which triggers circuitry 
to demultiplex that signal to the XA-1 analog output (giving a measure of the deposited 
energy) and to indicate which detector was struck. While this data is digitized by circuitry 
associated with the sandwich, further signals from the 128 detectors of that sandwich are 
ignored. This ‘deadtime’ has been measured at tdead = 98(3)/is. Thus the maximum com¬ 
bined rate at which the 256 sandwiches of BAT can detect photons is the saturation rate 
Rsat = ‘i^Q/tdead = 2.61(8) X lO^S-h 

If the BAT is measuring photons at a rate of Tmeas, then the average sandwich is dead 
for Tmeas/Rsat of the time, so the detector rate can be deadtime corrected by dividing by the 
livetime fraction: Tcorr = Tmeas/il - Tmeas/Rsat)- As Tmeas approaches Rsat the Uncertainty 
ill tdead makes T^orr poorly determined. Thus the deadtime correction from BAT measured 
count rate to photon rate is currently highly uncertain above Tcorr ~ ‘2-0Rsat- Above Tcorr ~ 
30i?sat the formal error range includes inhnite flux as the deadtime is consistent with 100%. 
The peak fluxes from SOPA indicate that the true photon rate on the BAT detector was 
~ 10^Rsat, and so even a much more accurate tdead determination would have resulted in a 
a poor estimate of Tcorr at the peak. Thus SOPA and ESP are needed to provide the peak 
measurements that BAT cannot. 
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Figure la: The SGR spike and tail light curve from the Burst Alert Telescope (BAT) on 
Swift at measured energy >50keV (64 ms bins). Although BAT was pointed 105° away 
from the SGR at the time of the main spike, it recorded 7 —rays above 60keV passing 
through and scattering within the spacecraft body and instrument shielding. As part of 
a pre-planned observing schedule, Swift slewed to observe a different source shortly after 
the main peak, reaching a steady pointing direction 61° from the SGR at 143 seconds. The 
spacecraft re-orientation improved the detection efficiency of the SGR, visible as an apparent 
(not intrinsic) rise in the light curve to a peak at 140 s. This is followed by a second slew to 
67°. 
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Figure lb: BAT deadtime-corrected count rate (all energies) during the complex leading 
edge of the main spike. (Note that horizontal scale is 10"^ larger than la.) Uncertainties in 
the deadtime correction (discussed in Supplementary Methods) make corrected count rates 
increasingly unreliable above 5 x counts per second. Time bins of lOOps are equivalent 
to the light-crossing-time of a neutron star diameter. More detailed lightcurves are shown 
in the supplementary figures section. 
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Figure 2: The pulse profile evolution of the magnetar SGR 1806-20 during the giant flare 
of 2004 December 27. Time through the flare increases from top to bottom. Each panel 
displays the pulse profile folded over four pulse cycles at one of four different time intervals 
during the flare. The times denoted at the top of each panel indicate the midpoint of the 
interval relative to the onset of the flare. Dnring the first half of the tail (~170 s), the peak 
centered at phase 0.7 grows in amplitnde as the primary peak fades until the two are nearly 
eqnal in height. Thereafter, the two peaks decay in lockstep while the relative amplitnde 
of the third peak at phase 0.2 increases. Overall, the pulse prohle becomes less sinusoidal 
during the course of the flare, that is, the power in the higher harmonics increases relative to 
the power at the fnndamental frequency, opposite to what was seen in SGR 1900+14 dnring 
the Angust 27 flare.The phases of the SGR 1806-20 pnlse peaks remain fixed, indicating 
a finalized magnetic geometry, and snpporting the notion that after the first spike no new 
magnetic energy is released, and only the trapped fireball energy leaks ont. The individnal 
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Supplementary Figure 1: BAT lightcurve of the precursor and tail of the giant flare with 
64 ms resolution. Successive plots are phase aligned at the 7.56 s neutron star spin frequency, 
with one cycle of overlap from plot to plot. The precursor, which occurred 142 s before the 
spike, is plotted at the appropriate spin phase. 
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plementary Figure 2: BAT deadtime-corrected lightcurve leading up to and during 
decline of the main spike, with 1 ms resolution. Corrected rates above 5 x 10^ have large 
irtainties in the deadtime correction. 
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Supplementary Figure 3: A spectral fit to SOPA data. This is an example of a spectral 
fit to data from one of the SOPA detectors to the peak 0.160 s of the main spike. The 
dashed line is the spectral model in photons, with a low energy power law a= -0.2 and an 
exponential cnt-off kT = 0.48 MeV. 







